General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



,1 


i ! O 


• 0 b-K ft £ £ 0 £ L l C « !! fi H I C i ft M P 


i. 1 k £. 


AEROSPACE SYSTEMS DIVISION 


mu ci e » m i n o ft i a k 


HOUSTON, TEXAS 7700ft 


TELEPHONE (AREA COOE 713) 4 ft ft • 0 0 ft 0 

Ref.: 642-2042 

Contract NAS 9-12200 
Job Order 73-763 


TECHNICAL MEMORANDUM 


ANALYSIS OF PRINCIPAL COMPONENT TRANSFORMED 

I 

LANDSAT DATA 


By 

R. A. Abotteen 


Approved By; 


"fJatto available under NASA sponsorship 
4n the interest of early and wide dts* 
semination of Earth Resources Survey 

information and without liability 
tor a fry we - made thereof.'-’ 

J. c. 

T. C. Minter , Supervisor 
Techniques Development 
Section 


Distribution: 

JSC/L. F. Childs 

J. D. Erickson (3) 

LEC/M. L. Bertrand 
B. L. Carroll 
J. E. Davis 
J. A. Quirein 
D. G. Saile 
P. C. Swanzy 
R. E. Tokerud 
Job Order File 
REDAF (3) 

Technical Library { 2 ) 


Original photography may be gircbtsed km 
tHUS Data Center 

Sioux Falls, SQ.^ 5 “71 yy . ± 


^ A 

DEG 1978 % 




\d> M 


|E7 9-10076.) ANALYSIS OF PRINCIPAL COMPONENT 
TRANSFORMED LANDSAT DATA (Lockheed 
Electronics Co.) 25 p HC A02/MF A01 

CSCL Q2C 

August 1976 


&HISIHAL GONTAiNS 

mvm mmmim 


N 79-15355 


63/43 


unclas 

00076 


LEC-9003 


S U B SID I A R Y OF LOCKHEED AIRCRAFT C 0 RP 0 RAT 



CONTENTS 


Section Page 

1. INTRODUCTION 1 

2. APPROACH 3 

2.1 IMPLEMENTATION OF THE PCOMP PROGRAM 3 

2.2 EXAMPLES : 3 

2.3 GENERATION OF PCOMP IMAGERY ........... 3 

3. THE EFFECT OF TRANSFORMATION ON CLASSIFICATION AND 

SIGNATURE EXTENSION PERFORMANCE 16 

4. CONCLUSIONS 20 

5. REFERENCES , 21 


i 



TABLES 

} 

Table Page 

1 EIGENVALUES, EIGENVECTORS, AND THEIR PROJECTION ONTO 

THE KAUTH SPACE OF SINGLE-PASS SITES 4 

2 EIGENVALUES, EIGENVECTORS, AND THEIR PROJECTION ONTO 

THE KAUTH SPACE OF TWO-PASS SITES 5 

3 CLASSIFICATION ACCURACIES FOR SEGMENT 1,183 WITH 

VARIOUS STATISTICS AND DATA TRANSFORMATIONS 17 

4 WHEAT AND NONWHEAT PERCENTAGES IN SEGMENT 1183 WITH 

VARIOUS STATISTICS AND DATA TRANSFORMATIONS 18 


ii 



FIGURES 


Figure Page 

1 Eight principal component images of the two-pass 
Dickinson site 

* < i 

(a) PCOMP 1 . . . 7 

(b) PCOMP 2 7 

(c) PCOMP 3 .......... 8 

(d) PCOMP 4 . . J . . ' 8 

(e) PCOMP 5 ............ . H 

(f) PCOMP 6 ....... 9 

(g) PCOMP 7 10 

(h) PCOMP 8 10 

2 First five principal component images of th« two- 
pass Morton site 

(a) PCOMP 1 11 

(b) PCOMP 2 . 11 

(c) PCOMP 3 12 

(d) PCOMP 4 .... • 12 

(e) PCOMP 5 13 

3 False color images of Dickinson site 

(a) Landsat date: 1454 14 

(b) Landsat date: 1634 . 14 

4 First three PCOMP false color image of Dickinson 

site (Landsat dates: 1454, 1634) 15 


iii 


1 . INTRODUCTION 


The Landsat multispectral scanner provides data in the form of 
four spectral images of the particular area of the Earth under 
observation. Ready and Wintz (ref. 1) have shown that the prin- 
cipal component (PCOMP) transformation applied to airborne- and 
satellite-gathered multispectral data is very useful for informa- 
tion extraction, since the first few PCOMP images contain essen- 
tially all the information present in the original spectral bands 
In this report, a PCOMP transformation was applied to single-pass 
and multipass Landsat data. The transformed data were analyzed 
visually and automatically. The PCOMP transformation is 

Y = MX (1) 

where 

X = vector of n spectral intensities associated with each pixel. 

M = n x n unitary matrix derived from the mixture covariance 

matrix Z x of the spectral bands such that the rows of M are 
the normalized eigenvectors of Z . 

X 

Y = vector of n principal components. 

The covariance matrix of the PCOMP transformed data then becomes 

Z = MZ M T 
Y X 



where X^, ••*, A n (the variances of the principal components) 

are the eigenvalues of Z v ordered such that A, > A,, > ••• > A 
ip ■ 12 n 

and M is the transpose of M. 
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An implementation plan for the program which outputs eigenvalues 
and eigenvectors is presented in section 2. Several examples and 
applications are also presented in section 2. Section 3 deals 
with the effect of the transformation on signature extension/ 
while section 4 states the conclusions 'drawn from this study. 
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2. APPROACH , 

• t 1. « 

The PCOMP transformation presented in this report was calculated 
as follows. A field containing all the picture elements (pixels) 
in a site was defined. An n-dimensional' covariance matrix was 
calculated for that field. The eigenvalues and normalized eigen- 
vectors corresponding to the covariance matrix were calcualted. 
The rows of the transformation matrix M were taken to be the 
normalized eigenvectors such that the first row corresponds to 
the largest eigenvalue and the second row corresponds to the 
second largest eigenvalue, etc. 

2.1 IMPLEMENTATION OF THE PCOMP PROGRAM 

The University of Houston SYMAT program (ref. 2) was modified 
and implemented on the Laboratory for Applications of Remote 
Sensing (LARS) terminal under the file name PCOMP. The modifi- 
cations made allowed the PCOMP program to read the covariance 
matrix from LARSYS-punched cards. It also added the capability 
of projecting the principal components onto the Kauth (ref. 3) 
space. For multipass data, the principal components were pro- 
jected onto the Kauth space four channels at a time since the 
Kauth transformation is four-dimensional. 

2.2 EXAMPLES 

The eigenvalues and eigenvectors of for some U.S. Department 
of Agriculture (USDA) Statistical Reporting Service (SRS) and 
intensive test sites were computed using single-pass and multi- 
pass data. The sites used were Marion, Dickinson, and Morton 
counties in Kansas. The' results are given in tables 1 and 2. 

Each table also shows the projection of the eigenvectors onto 
the Kauth space . The first two principal components in table 1 
projected onto the Kauth space are primarily Bright (first com- 
ponent) stuff and Green (second component) stuff. The remaining 
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TABLE 1.— EIGENVALUES, EIGENVECTORS, AND THEIR PROJECTION ONTO 
THE KAUTH SPACE OF SINGLE-PASS SITES 


Eigenvalues 

of 


34.4 

13.6 

1.4 

0.8 


54.2 

49.8 

4.3 

0.9 



77.5 

52.4 

3.1 

1.0 


Eigenvectors of 

E 

X 

0.42 

0.56 

0.65 

0.54 

0.30 

i 

o 

• 

u> 

-0.63 

0.44 

0.78 

-0.54 

0.10 

-0.29 

-0.32 

0.06 

0.52 

-0.79 

-0.03 

0.12 

-0. 83 
0.13 


0.34 

0.93 

-0.01 

-0.91 

0.34 

-0.07 

0.24 

• 

O 

1 

0.01 

0.53 

-0. 81 

isa 

0.50 

0.70 

0.37 

MM 

-0.70 

0.45 

0.44 

-0.75 

0.50 

-0.20 

0.40 

0.45 

-0.11 

-0.52 

0.72 

mm 

-0.29 

0.78 

0.55 


-0. 83 

-0.31 

-0.06 

00 

• 

0 

1 

0.47 

-0.09 

0.24 

00 

(N 

• 

0 

1 

0.01 

0.54 

-0.80 


Projection of eigenvectors 
onto Kauth space 


0.35 

-0*38 

• 0.11 

0.85 


0.68 

-0.52 

0.23 

-0.47 


0.99 

• 0.11 

•0.03 

■ 0.01 


0.98 

- 0.21 

- 0.02 

0.01 


0.38 

■0.92 

■0.06 

■ 0.01 


-0.44 

-0.34 

0.82 

0.14 


I 


.18 

.98 

.01 

0.00 


0.11 

0.99 

• 0.01 

• 0.02 



0.22 

0.98 

0.01 

-0.03 


0.92 

0.38 

0.04 

0.01 


■0.98 

-0.18 

• 0.01 

■ 0.02 


- 0.02 

0.02 

-0.99 

0.12 


0.01 

0.20 

0.98 

0.04 


0 
6 

1.00 

0.06 


0.00 

0.00 

0.39 

■0.92 


0.00 

0.01 

0.06 

- 1.00 


.01 
.01 
.92 
0.39 








































Kauth- projected PCOMP eigenvectors are primarily Yellow (third 
component) stuff and Nonesuch (fourth component) or a combination 
of them. In table 2, the first four Kauth-projected PCOMP eigen- 
vectors are primarily Bright and Green stuff while the remaining 
projections are primarily Yellow stuff and Nonesuch. 

2.3 GENERATION OF PCOMP IMAGERY 

Landsat data were transformed using the unitary matrix M defined 
in equation (1). The transformation was performed using the 
LARSYS data transformation processor. The processor transformed 
then rescaled the data to lie between 0 and 255. Black-and-white 
images of principal components 1, 2 , 2 , 4, 5, 6, 7, and 8 corre- 
sponding to the two-pass SRS Dickinson site (Landsat dates: 1454, 

1634) are shown in figure 1. Figure 2 shews the first five prin- 
cipal component images of the Morton intensive test site (Landsat 
dates: 1454, 1727). There is essentially no information in 

PCOMP 5 or higher. This indicates that the inf ormation:: in the 
Landsat bands can be packed in fewer bands using the principal 
components that correspond to the largest eigenvalues of I . 

A false color image of the first three principal components was 
generated using multitemporal data. Figure 3a shows a false 
color image of the original Landsat bands of Dickinson in the 
first pass (Landsat date: 1454)* Figure 3b shows a false color 

image of Dickinson using another pass (Landsat date: 1634) . 

The first three principal components of the eight-channel 
Dickinson site (Landsat dates: 1454, 1634) were used to generate 

the false color image shown in figure 4. 
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(a) PCOMP 1. 


(b) PCOMP 2. 

Figure 1. — Eight principal component images of the 
two-pass Dickinson site. 
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Figure 1. — Continued. 
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(e) PCOMP 5. 
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(f) PCOMP 6. 


Figure 1. - Continued. 





(h) PCOMP 8 
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(a) PCOMP 1. 



(b) PCOMP 2. 


Figure 2. — First five principal component images of the 

two-pass Morton site. 
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(c) PCOMP 3 


(d) PCOMP 4 


Figure 2. — Continued 
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(e) PCOMP 5. 


Figure 2. — Concluded. 
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Figure 3. - False color images of Dickinson site. 
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Figure 4. — First three PCOMP false color image of 
Dickinson site (Landsat dates: 1454, 1634). 
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3. THE EFFECT OF TRANSFORMATION ON CLASSIFICATION AND 
SIGNATURE EXTENSION PERFORMANCE 

1 >! 

A consecutive-day acquisition (CDA) for Large Area Crop Inventory 
Experiment (LACIE) segment 1183 was chosen to study the effect 
of PCOMP transformation on classification and signature extension 
performance. This segment was acquired in the crop year 1975. 

The image with acquisition date 75110 was clear, while the 75111 
acquired image was hazy. This segment along with several other 
CDA segments was used previously for testing signature extension 
algorithms (refs. 4 and 5). This segment was chosen because 
there is a large difference in haze level between the two 
acquisitions . 

As mentioned previously in this report, the modified LARSYS data 
transformation processor transformed the data then rescaled it 
to have a range from 0 to 255. Various combinations of data 
manipulations were applied to both acquisitions of segment 1183. 

The PCOMP transformation used for both acquisitions of segment 1183 
was derived as follows. Both acquisitions were merged to form an 
image of 234 lines with 196 pixels in each line. Each acquisition 
has four channels, and thus the merged image also has four chan- 
nels. A mixture covariance matrix Z c was calculated using all 
the pixels in the merged image. Eigenvalues and eigenvectors of 

Z were calculated using the PCOMP program. The PCOMP transfor- 
c 

mation was the matrix of the first two eigenvectors of Z . Once 
the data in each acquisition of segment 1183 were transformed, 
calculations of statistics and all classifications were performed 
in the two-dimensional transformed space. Each acquisition of 
segment 1183 was classified using local statistics and statistics 
from the other acquisition. The classification results are shown 
in tables 3 and 4. These tables also show some results obtained 
from reference 5 in which the original Landsat channels used for 
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75111 75110 Acquisition 


TABLE 3.- CLASSIFICATION ACCURACIES FOR SEGMENT 1183 WITH VARIOUS STATISTICS 

AND DATA TRANSFORMATIONS 


Field accuracy, % 


Crop 


Original Landsat 
data 


ROOSTER 

R(S) 



Wheat 89 . 1 

Non- 95 . 6 

wheat 

Wheat — 

Non- — 

wheat 



Transformed Landsat data 


UT 

PCOMP/SC 

PCOMP/SC 

SC 

SC 

PCOMP/no 

PCOMP/no 

local 

UT 

local 

UT 

SC local 

SC UT 

6.6 

&8. 3 

86.9 

89.1 

68.4 

89.2 

4.5 

44.8 

96.0 

94.8 

94.8 

94.4 

95.2 

64.3 

— 

88.2 

88.3 

89.5 

53.6 

88.6 

10.7 

— 

93.3 

89.7 

91.3 

58.7 

94.4 

80.1 


NOTES: ROOSTER R(S) = Rank Order Optimal Signature Transformation Estimation Routine- 

that Uses subclass means. 

UT =" untransformed statistics from training segment used in classification. 
PCOMP/SC = principal component transformed and rescaled data. 

SC = rescaled data. 

Local = local statistics used in classification. 
























75111 75110 Acquisition 


TABLE 4.- WHEAT AND NONWHEAT PERCENTAGES IN SEGMENT 1183 WITH VARIOUS STATISTICS 

AND DATA TRANSFORMATIONS 






Proportion, 

% 




Crop 

Original Landsat 
data 


Transformed Landsat 

data 



Local 

ROOSTER 
R (S) 

UT 

PCOMP/SC 

local 

PCOMP/SC 

UT 

SC 

local 

SC 

UT 

PCOMP/no 
SC local 

PCOMP/no 
SC UT 

Wheat 

55.1 

43.1 

9.1 

53.7 

49.2 

54.7 

30.8 

54.0 

8.6 

Non- 

wheat 

35.2 

39.6 

22.6 

37.8 

42.5 

35.8 . 

55.9 

39.1 

26.9 

Wheat 

— 

_ 

— 

51.2 

55.7 

52.0 

40.1 

50.4 

4.3 

Non- 

wheat 

— 

— 

_ 

40.7 

36.0 

39.4 

11.9 

41.1 

67.7 


NOTES: ROOSTER R(S) = Rank Order Optimal Signature Transformation Estimation Routine 

that uses subclass means. 

UT = untransformed statistics from training segment used in classification 
PCOMP/SC = principal component transformed and rescaled data. 

SC = rescaled data. 

Local - local statistics used in classification. 




classification purposes and the training statistics were trans- 
formed with the Rank Order Optimal Signature Transformation 
Estimation Routine (ROOSTER) . 

Table 3 shows the training field accuracy for segment 1183 using 
original and transformed Landsat data. The local training field 
classification accuracy obtained by using the two-dimensional 
PCOMP transformed Landsat data is almost the same as the accuracy 
obtained by using the original Landsat data. The same observa- 
tion can be noted in table 4 when proportions obtained by using 
the two-dimensional PCOMP transformed and original Landsat data 
are compared. Based on this very good classification performance 
of PCOMP transformed data, the transformation can be used as a 
feature selector. 

The classification performances presented in tables 3 and 4 
indicate that PCOMP transformation followed by rescaling of the 
data resulted in excellent signature extension performance. 

When the PCOMP transformation was used alone, the signature 
extension performance was bad; however, the signature extension 
performance was fair when scaling of the data alone was used. 
Apparently scaling of the data corrected for some of the atmos- 
pheric differences between acquisitions 75110 and 75111, but the 
PCOMP/SC transformation corrected for all differences between 
the acquisitions. 
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4. CONCLUSIONS 

1. The first few principal components corresponding to the 
largest eigenvalues contain most of the information of a 
multispectral Landsat image. The local classification results 

indicate that the PCOMP transformation is an excellent feature 

1 

selector which is in agreement with the results obtained by 
Ready and Wintz (ref. 1) . 

2. Based on the examination of two analyst interpreters (AI's), 
the multi temporal principal component false color imagery 
may be helpful for selecting training fields in LACIE. 

3. A PCOMP transformation followed by a rescaling of the data 
gave excellent signature extension results; therefore, it is 
assumed that most of the atmospheric differences between 
acquisitions 75110 and 75111 of segment 1183 were eliminated. 
Apparently a significant part of this was due to scaling 
alone, since scaling of the original data caused a significant 
improvement in the PCOMP/SC signature extension performance. 
Further studies in this area are recommended. 
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